Accepted by ApJ on July 1, 2009 

Preprint typeset using 1^1^^ style cmulatcapj v. 10/09/06 



NEAR-INFRARED SPECTROSCOPY OF SDSS J0303-0019: 
A LOW LUMINOSITY, HIGH EDDINGTON RATIO QUASAR AT Z 

.t2.3 



6 



J. D. Kurk\ F. Walter' 



O 
O 

(N 



(N 

o 
u 

6 



> 
m 
p 

o 
o 



X 



X. Fan^'', L. Jiang^ S. Jester\ H.-W. Rix\ D. A. Riechers^'^ 

Accepted by ApJ on July 1, 2009 



ABSTRACT 

We present sensitive near-infrared VLT ISAAC spectroscopic observations of the z — 6.08 quasar 
SDSS J030331.40-001912.9. This QSO is more than a magnitude fainter than other QSOs at z - 6 for 
which NIR spectroscopy has been obtained to date and is therefore presumably more representative 
of the QSO population at the end of Cosmic Reionization. Combining rest-frame UV continuum 
luminosity with the width measurements of the Mg II and C IV lines, we derive a black hole mass of 
2io'5 X 10® Mq, the lowest mass observed for z ~ 6 QSOs to date, and derive an Eddington ratio of 
l-6lo'g, amongst the highest value derived for QSOs at any redshift. The Spitzer 24 //m non-detection 
of this QSO does not leave space for a significant hot dust component in its optical/near-infrared SED, 
in common with one other faint QSO at z = 6, but in contrast to more than twenty more z = 6 QSOs 
and all known lower redshift QSOs with sufficiently deep multi-wavelength photometry. We conclude 
that we have found evidence for differences in the intrinsic properties of at least one z ^ 6 QSO as 
compared to the lower-redshift population. 

Subject headings: Galaxies: high-redshift — early Universe — quasars: emission lines — quasars: 
individual: SDSS J030331.40-001912.9 



1. INTRODUCTION 

The radiation from luminous quasi stellar objects 
(QSOs) can currently be studied up to redshift 6.5 and 
thus allows one to probe, over much of the Universe's 
age, the properties of their central black holes (BHs, e.g., 
Willott ct al. 2003), their broad line region (BLR, e.g., 
Freudling et al.ll20p3[). and the ioni zation state of inter- 
galactic gas fe.g.. iFan et al.|[2006bD . The most distant 
QSOs are of importance for both the study of galaxy 
evolution and large-scale structure. Among the most 
surprising discoveries in this field are the detection of BH 
mass es > 10^ Mq (e.g.. iBarth et al.ll2003t IWillott et all 
l2003l ) and the apparent constancy of the observed QSO 
BLR properties with reds hift, even up to z — 6.4 
(e-g.- lPentericci et al.ir2002[ ). corresponding to an age of 
870 Myr. E.g the Fe II/Mg II ([a cudling ct al. 2003; 
Maiqlino etall 120031 : llwamuro etH. 2004: Kurk ct al.; 
20071 hereafter K07) and N V/C IV (Jiang et al. 2007.) 
emission line ratios observed in the BLR of z > 6 QSOs 
are as high as that observed in low-redshift objects. 
Thus, near-solar metallicity in the centers of some of the 
earliest systems and accretion of > 10^ solar masses in 
a central BH can be attained in less than 1 Gyr. One 
caveat, however, has to be taken into account: until now 
only the most luminous objects at z ~ 6 have been stud- 
ied. These objects may experience faster evolution than 
less luminous objects and therefore resemble low to inter- 
mediate redshift QSOs already at z ~ 6. Secondly, the 
relationships established using low-luminosity AGN (at 
low redshift) may not be valid for objects at the upper 
end of the luminosity distribution. As an example, the 
determination of the BH mass in active galaxies was cali- 
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brated with reverberation mapping observations of AGN 
that are one to two orders of mag nitude fainter than the 
z 6 QSOs observed until now (|Kaspi et al.ll2000f ). It 
is therefore important to extend z > 6 studies to lower- 
luminosity objects. Here we present the first observations 
of a QSO at z > 6 with luminosity lower than Lboi — 10"*^ 
ergs-i: the QSO SDSS J030331. 40-001912.9 (hereafter 
J0303-0019). This QSO was detected in the Sloan Digi- 
tal Sky Survey (SDSS) Deep Stripe with a magnitude of 
ZAB = 20.9. Its discovery spectrum shows a very strong 
and n arrow Lya emission line at z = 6.070 (jJiang et al.l 
l2008f ). As this line falls in the z-band, the actual con- 
tinuum luminosity is low, mi45o = 21.3, 1.5 magnitudes 
fainter than the mean of the z 6 QSO samples studied 
in K07 and iJiang et al.l ()2007t ) . In addition, this is one 
of the two out of about twenty z 6 QSOs that are not 
detected at 24 /xm with MIPS /Spitzer, implying a rest 
frame 3.5Atni /4400A ratio far lower than for any other 
known QSO (jJiang et al.ll200d L. Jiang et al. in prep.). 

2. OBSERVATIONS AND REDUCTION 

2.1. Near-infrared spectroscopy 

NIR spectroscopy of J0303-0019 in the Y and K 
bands was obtained in service mode on December 24 and 
November 26 and 28, 2007, respectively, with ISAAC 
(jMoorwood et al.lll998[) at VLT UTl. We used the short 
wavelength arm, equipped with a 1024x1024 Hawaii 
Rockwell array which has a pixel scale of 0.147". The 1" 
slit results in a spectral resolution of 550 and 450 for the 
Y and K bands, respectively. Observations were carried 
out in a standard ABBA offset pattern, with additional 
dithering to reduce the influence of bad pixels. Individ- 
ual exposures were 148s, repeated 16 times per observing 
block. Nine such blocks were obtained in K band, for a 
total exposure time of ~5.9hr, and two in Y band, for a 
total time of ~ 1.3 hr. The data were reduced using ESO's 
ISAAC pipeline software, which does flat fielding, back- 
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TABLE 1 

NIR Vega magnitudes of J0303-0019 



V 


J 


H 




20.60±0.14 


20.44±0.08 


19.78±0.08 


18.95±0.09 



ground subtraction, rectification and combination of the 
individual frames in one observing block. We used sky 
lines to determine the wavelength calibration, obtaining 
an absolute la error of 3 and 7 A, respectively, in Y and 
K band. One-dimensional spectra were extracted from 
the pipeline products in IRAF. B and G stars with known 
NIR photometry were used to correct for telluric absorp- 
tion and initial flux calibration of the one-dimensional 
spectra. Finally, the one-dimensional spectra were aver- 
aged to obtain the final spectra. These were flux cali- 
brated using the broad band magnitudes obtained with 
the Calar Alto 3.5m telescope (see Sec. l2.2|) . The reduc- 
tion is similar to that carried out by K07, where more 
details are reported. 

2.2. Near- infrared broad hand imaging 

For accurate flux calibration, the QSO SDSS J0303- 
0019 was i maged with the pr ime-focus NIR camera 
Omega2000 (|Kovacs et all I2004D at the 3.5m telescope 
at Calar Alto, Spain. This instrument has a 2kx2k 
IIAWAII-2 detector with a pixel scale of 0.45", result- 
ing in a field of view of 15.4' x 15.4'. J0303-0019 was 
observed on September 10 and 11, 2008, for 20, 20, 25, 
and 77 minutes in Y, J, iJ, and Kg bands, respectively. 
After dark subtraction and flat fielding with twilight ex- 
posures, the images were sky subtracted, registered and 
averaged with XDIMSUM^ in IRAF'^. The seeing full 
width at half maximum (FWHM) on the combined im- 
ages is 1.3", 1.4", 1.3", and 1.0", respectively. The com- 
bined images were astrometrized using the public ser- 
vice provided by astrometry.net ()Hogg et al.ll2008f ) us- 
ing a second order polynomial to correct for distortions 
(which are small for Omega2000). Subsequently, we used 
about 80 stars from the 2MAS S All-Sky Catalog of Point 
Sources (jSkrutskie et al.ll2006D to obtain the magnitude 
zero points. For F-band, we use d the UKIDSS second 
data release (jWarren et al.l [20071 ) catalog. Only half of 
the observed field is covered by UKIDSS, but this still 
provides 31 stars usable for the fiux calibration. Sys- 
tematic errors in the fiux calibration are estimated to be 
below 0.1 mag. The resulting fluxes (us ing MAG_BEST 
in SExtractor, iBertin fc Arnout"3 Il996l ) are reported in 
Table □ 

3. RESULTS 

3.1. Optical and NIR spectra 

In the top panel of Fig. [U we show a clos e-up version 
of the optical spectrum, flrst presented in iJiang et al.l 
(|2008D . taken at the W.M. Keck Observatory. For com- 
pa rison, the SDSS QS O composite spectrum published 
bv iVanden Berk et al.l (2001, dashed curve), normalized 
to emission line free regions of our observed spectra, is 
overlayed, showing that the emission lines of J0303-0019 

6 XDIMSUM is a variant of the DIMSUM package developed by 
P. Eisenhardt, M. Dickinson, S.A. Stanford, and J. Ward. 

IRAF is distributed by NOAO, which are operated by AURA, 
Inc., under cooperative agreement with NSF. 
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Fig. 1.— Top: Optical spectrum of J0303-0019 (black his- 
togram), the same spectrum smoothed over 20 A (thick grey curve) 
and the noise per pixel (blue line). A simultaneous fit (thick green 
curve) of the N V line and the broad and narrow components 
of the Lya line is over plotted on top of a power-law with slope 
of ax = —1.5, normalized to a line free region of the spectrum 
(short-dashed line). In addition, the wavelength of the onset of 
the GP trough and of Lya at the systemic redshift are indicated 
by vertical dashed lines. The normal ized SDSS QSO composite 
spectrum by IVanden Berk et al.l l|2001l ) is over plotted in red for a 
visual comparison of the line widths. Middle and bottom: NIR Y 
and K band spectra of J0303-0019, including the C IV and Mg II 
lines, respectively, smoothed over 5 pixels. The grey regions are 
seriously affected by sky lines and omitted from the fit. Lines and 
colors as in top panel. The slope of the power-law continuum here 
is ax = —1.0. 

are much stronger and narrower than those present in 
the composite spectrum. We have made a simultane- 
ous fit of the broad and narrow Lya and N V lines 
in the optical spectrum, after subtracting a standard 
(IVanden Berk et al.|[200l power law co ntinuum with a 
slope Fa cx A'^-^ as in lJiang et a"n(|2008D . We fit only the 
spectrum redwards of the strong absorption feature seen 
in the Lya line. The best fit parameters are listed in Ta- 
ble m and the best fit curves are over plotted in Fig. [TJ 
The equivalent width (EW) and FWHM of N V mea- 
su red by us ar e abou t 20% higher than those reported 
bv I Jiang et al.l ()2006( ). who fixed the central wavelength 
of all fitted lines to the same redshift. 

In the middle and bottom panels of Fig. [l] we show 
the VLT ISAAC Y- and iiT-band spectra, containing the 
C IV and Mg II emission lines, respectively. There is also 
a tentative detection of Si IV/0 IV] 1398 at 9900 A in the 
V-band, which we do not try to fit. The obtained NIR 
photometry of J0303-0019 (Tabled]) is consistent with a 
power law SED with slope a a = —1.5, as measured in 
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TABLE 2 

J0303-0019 Emission line measurements 



Line 


a 


z 


-'^central 


FWHM'' 


Flux" 


EWo 


Mbh 










A 


A 


kms ^ 


10-1^ cgs 


A 


10** Mq 




Ly« (broad) 


-1.5 


6.073±0.001 


8599±1 


125± 5 


4350±165 


255± 6 


87±2 






Lya (narrow)"^ 








29± 2 


lOOOit 70 


109±11 


37±4 






N V 




6.063±0.001 


8758±1 


81± 2 


2760± 80 


90± 2 


32±1 






ClV 


-1.5 


6.071±0.001 


10953±1 


65± 2 


1780± 65 


54± 2 


37±1 


1.7±0.1 


1.8±0.2 


CIV 


-1.0 


6.071±0.001 


10953±1 


65± 3 


1772± 70 


54± 2 


36±1 


1.6±0.1 


1.9±0.2 


CIV 


0.0 


6.071±0.001 


10953±1 


64± 3 


1760± 70 


53± 2 


35±1 


1.5±0.1 


2.0±0.2 


Mgii 


-1.5 


6.078±0.001 


19804±4 


149± 9 


2261±130 


16± 1 


25±1 


2.3±0.3 


1.3±0.2 


Mg II 


-1.0 


6.078±0.001 


19803±4 


154± 8 


2340±125 


17± 1 


27±1 


2.5±0.3 


1.2±0.1 


Mg II 


0.0 


6.077±0.001 


19801±3 


167± 8 


2525±120 


19± 1 


32±1 


3.0±0.3 


l.OitO.l 



lower redshift QSOs (|Vanden Berk et"aLll200ll ). Imper- 
fect flux calibration of the NIR spectra, however, may 
affect the slope of the reduced Y and K band spectra. 
Indeed, the slope observed in the K band seems redder 
and may be fit better by a slope of = —1.0 or even 
a\ = 0. As the NIR spectra are too short and the data 
not of sufficient quality to constrain the slope, we carry 
out emission line fits employing a range of slopes a\ fixed 
to = 0, -1.0, and -1.5, to cover the complete range of plau- 
sible slopes. We take the slope of a a — —1.0 as the best 
compromise and use the parameters derived from fits em- 
ploying the two extreme values to construct conservative 
error estimates on the fitted parameters. The C IV and 
Mg II lines are fit by single Gaussian functions on top 
of the co-addition of the power-law continuum, a Balmer 
pseudo-continuum and a template of Fe II lines (see also 
iKurk et aTll2007t l. The strength of the latter two com- 
ponents were scaled with the power-law continuum, ac- 
cording t o their proportions i n the composite QSO spec- 
trum by IVanden Berk et al.l ()2001l ). As the strength of 
the Fe II emission is not well constrained by the current 
spectra, but clearly lower than four times the strength 
in the SDSS composite spectrum, we have also carried 
out the fitting assuming a four times stronger and four 
times weaker Fe II contribution to establish proper error 
estimates. 

The resulting best-fit parameters are displayed in Ta- 
ble H The C IV and Mg II lines have (rest frame) 
FWHMo less than 2000 and 2600 kms-\ and EWq of 
36±lA and 27t2 A, respectively. The mea n FWHMp 
of the se lines in the SDSS QSO sample of IShen et al] 
(|2008l hereafter SOS) is 5700 and 5850 kms-\ respec- 
tively. Only 2% and 1% of the SDSS QSOs have C IV 
and Mg II FWHMq smaller than 2000 and 2600 kms^^, 
respectively. 

3.2. Black hole mass measurements 

The widths of the Mg II and C IV lines, combined 
with an independent estimate of their characteristic dis- 
tance from the center, can be used as a measure of the 
mass of the black hole (Mbh)- The relation between line 
width, luminosity and Mbh has been calibrated empiri- 
cally at low redshift using reverber a.tion mapping (e.g ., 
iKaspi et al.ll20000 . As in K07 and iJiang et all (|2007h . 
we use the relations for C IV and Mg II presented by 
IVestergaard fc Petersoiil (|2006f ) and iMcLure fc Dunlod 
( 2004[ ). respectivelv. which are also those applied to a 
sample of more than 60,000 SDSS QSOs by SOS. Due to 
the scatter of individual data points around this relation 
Mbh can be determined only within a factor of three. 



For J0303-0019, we derive Mbh with these methods and 
take into account the range of power-law slopes and Fe II 
template strength described in Sec. 13.11 We obtain Mbh 
of l.Qt^^l and 2.5tg|xl0S Mq, respectively, from the 
C IV and Mg II lines. These are the lowest Mbh derived 
in a 2; > 6 QSO to date, showing that, with this sample of 
fainter z ~ 6 QSOs, we are indeed probing down the mass 
function of BHs. We note that the central wavelength of 
the C IV line is blue shifted by 350 kms^^ w.r.t. the 
Mg II fine. According to SOS, this modest shift (< 1000 
kms^^) implies that the C IV line profile is dominated 
by gravitational effects and that Mbh can therefore be 
reliably estimated from this line. 

The Mbh derived from the line widths and correspond- 
ing Eddington luminosities of logLEdd = 46.3lto;o9, 
and 46.511qo7 ergs~^, respectively, are small com- 
pared to the QSO's bolometric luminosity of logLboi = 
46.59±0.04 ergss^^, implying a high Eddingtion ratio 
(rEdd)- From optical, IR, and radio data (L. Jiang et al, 
in prep.), the bolometric luminos ity of J0303-0019 is de- 
termined as in iJiang et al.l ()2006f ) to be logLboi — 46.59 
ergss"^. It is close to the values of logLboi = 46.62, 
and 46.65 ergs~^ obtained using the monochromatic lu- 
minosities at 1350 A and 3000 A, respectively, and the 
bolometric correction factors^ employed by SOS (BC1350 
= 3. SI, BC3000 = 5.15). Reconciling these values, we ob- 
tain the estimated bolometric luminosity stated above. 
The Eddington ratio of most QSOs is between 0.01 and 
1 and only S04 out of the 62,1S5 SDSS with Mbh deter- 
mined from their line-width have rEdd> 1 (SOS). The 
Eddington ratios for J0303-0019 derived from the C IV 
and Mg II lines are rEdd = l-9to4 and I.2I02 (see Ta- 
ble [2] and Fig. [2]). Less than 2% and 0.2% of the SDSS 
QSOs have such an rEdd (C IV) > 1.9 and rEdd (Mg II)>1.2. 
Even if we assume that the other four (faint) z ~ 6 QSOs 
of the compl ete sample of QSO s discovered in the SDSS 
Deep Stripe (jJiang et al.ll2008D have rEdd <1.2, the frac- 
tion of Lboi < 10*^ ergs~^ QSOs at z '--^ 6 with rgdd >1.2 
is 33%, almost 40 times higher than at lower redshifts. 

3.3. Abundance indicators 

The chemical composition of gas clouds in the BLR 
can be explored by comparing observed emission line 
ratios to models. The heavy elements in the BLR are 
products of the star formation history in the host galaxy 
and therefore indicators of p ast star formation in at least 
the very center of the host. iNagao et al.l (|2006f l present 

* These bolometric correction factors are multiplicative in flux, 
unlike the common definition in units of magnitude. 
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Fig. 2.— The accretion luminosity of J0303-0019 and J0005- 
0006: the figure shows the Mbh ^ Lbol diagram reproduced from 
SOS, including ~60,000 SDSS QSOs at 0.1 < ^ < 4.5. Colors 
show Mbh estimates using the following virial estimators: red for 
H/3, green for Mgll, blue for CIV. Three diagonal lines show 0.01, 
0.1, and 1 Lfidd, as indicated. Superposed are the data points for 
J0303-0019 (this work) and for J0005-0006 (K07). Here, we have 
applied the cosmology specified in SOS, resulting in Mbh (Lbol) 
larger by 5% (7%), but r^dd similar to those in our cosmology. 

several emission line ratios observed in SDSS QSOs at 
2.0 > z > 4.5 and m odels explaining the ratios. Accord- 
ing to iNagao et al.l . the strength of the N V/C IV ratio 
is strongly correlated with the luminosity of the QSO, 
which can be interpreted in terms of higher gas metal- 
licity in more luminous QSOs. Although J0303-0019 is 
one of the faintest QSOs at z ~ 6, it has a Mb of -26.2, 
which makes it moderately luminous in c omparison with 
the SDSS QSOs studied bv iNagao et"all . We measure a 
ratio of N V/C IV — 1.61±0.07, the highest ratio among 
the known z ~ 6 QSOs, and also higher than those mea- 
sured in the composite spectra of quasars in the luminos - 
ity range -25.5 > Mb > -26.5 by Nag ao et al.l (|2006D . 
This suggests that the BLR has a high meta l abundance, 
possi bly as high as IOZ0 (see Fig. 29 in INagao et alJ 
|2006[) . but the authors warn that estimates of the BLR 
metallicity using only emission line ratios involving N V 
might be quite uncertain (see also Ham ann et al.l[200^ . 

3.4. Stromgren sphere 

The intense UV radiation of the QSO keeps a region 
around the QSO ionized, the Stromgren sphere, which 
will actually only be spherical if the IGM surrounding 
the QSO is homogeneous. Outside this region, UV pho- 
tons are absorbed by neutral hydrogen and this causes 
complete absorption blueward of the wavelength of red- 
shifted Lya, the so-called GP trough (Gunn fc Pe tersonl 
[1^5). For reasons discussed in Fan et a l. (2006_d), the 
Stromgren sphere radius is difficult to measure and we 
rather use the size of the proximity zone (Rp), defined 
as the region in which the transmitted ffux is >10% of 
the (extrapolated) continuum emission when smoothed 
to a resolution of 20 A (as shown by the grey line in 
Fig. [1]). We measure the blue edge of this region to be 



at 8513±1 A (dashed line in Fig.[T|). Combined with the 
systemic redshift of z — 6.078 derived from the Mg II 
line (redmost dashed line in Fig. [T]), we compute (K07) 
the co-moving radius of the Stromgren sphere to be Rs 
= 4.5±0.3Mpc. This is a ty pical va l ue for luminous 
z ~ 6 QSOs (see, e.g., K07. iWalter et al.il2003HFan et all 
l2006al) . The average size of QSO Stromgren spheres, 
however, decreases by a factor ~2.5 over 5.7 < z < 6.4, 
consiste nt with an increas e of the neutral fraction with 
redshift (|Fan et al.ll2006bD . 

4. DISCUSSION AND CONCLUSION 

SDSS J0303-0019 has the lowest BH mass known at 
z 6 and is a truly peculiar QSO: it has one of the 
highest Eddington ratios observed for a QSO, it has the 
highest N V/C IV ratio observed in z ^ 6 QSOs, which 
is also significantly higher than that observed in lower 
redshift composite spectra of QSOs, and it is one of the 
two z '--^ 6 QSOs that is not detected in deep 24 //m ob- 
servations (jJiang et al.l l2006l . L. Jiang et al. in prep.). 
The other QSO which is not detected at 24 ^m is SDSS 
J000552.34-000655.8 (hereafter J0005-0006), which also 
has a high Eddington ratio: 2.4 as derived from the 
Mg II based Mbh and 1.0 as derived from the C IV based 
Mbh(K07). Detections at 250 GHz and 1.4 GHz have also 
been attempted but neither of the objects were detected 
down to rms values of 0.51, 0.48 mJy and 62, 130 /iJy, 
respectively (Wang et al. 2008,). In contrast, these QSOs 
were detected at 3.6, 4.5 and 5.8/xm with IRAC/Spitzer 
piang et al.l I2006I . L. Jiang et al. in prep.), implying a 
rest frame 3.5/xm/4400A ratio far lower than for any ob- 
ject in the z ~ 6 QSO sample, and lowe r than f or any 
of t he PC QSOs observed by [Neugebauer et all ()1987D 
and lHaas et al.l (|2000l [20031) for which the sensitiv ity at 
3.5mu rest-frame is not an issue. I Jiang et all ()2006f ) show 
that the optical/IR SEDs of the observed z ~ 6 QSOs 
can be fit by two components: a power-law representing 
an accretion disk and a black-body representing hot dust. 
The respective (Icr) upper limits of 21/xJy and 9/iJy at 
24 /xm for J0303-0019 and J0005-0006, significantly lower 
than the expected flux of ~160/iJy based on SEDs of 
other z '-^ 6 QSOs, suggests the absence of hot dust. 
The implied lack of dust could be explained if the gas 
is metal-poor, as in a galaxy host which has not experi- 
enced major star formation yet. However, this interpre- 
tation would contradict the potentially high metallicity 
based on the measured N V/C IV ratio. Alternatively, 
the non-detection of the hot dust component could be 
due to different dust properties, or due to a surrounding, 
cool dust torus hiding the hot component (jJiang et al.l 
I2OO6I) . 

We conclude that, with the observation of QSOs with 
bolometric luminosity < lO^'' ergs~^, we start to probe a 
population of QSOs that shows different observed prop- 
erties as compared to higher luminosity z ~ 6 QSOs and 
lower-redshift QSOs of all luminosities. We can only 
speculate that the hosts of the lower-luminosity QSOs 
at z ~ 6 may be younger than higher-luminosity QSOs 
at the same redshift because the host galaxies of the lat- 
ter are more massive and have therefore formed earlier. 
It is, however, unclear at this stage what physical mech- 
anisms are causing the observed differences. 
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